Nonlinear Band Gap Transmission in Optical Waveguide Arrays 
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The effect of nonlinear transmission in coupled optical waveguide arrays is theoretically investi- 
gated via numerical simulations on the corresponding model equations. The realistic experimental 
setup is suggested injecting the beam in a single boundary waveguide, linear refractive index of 
which (no) is larger than one (n) of other identical waveguides in the array. Particularly, the effect 
holds if Lu(no — n)/c > 2Q, where Q is a linear coupling constant between array waveguides, ui is a 
carrier wave frequency and c is a light velocity. Making numerical experiments in case of discrete 
nonlinear Schrodinger equation it is shown that the energy transfers from the boundary waveguide 
to the waveguide array above certain threshold intensity of the injected beam. This effect is ex- 
plained by means of the creation and propagation of gap solitons in full analogy with the similar 
phenomenon of nonlinear supratransmission [F. Geniet, J. Leon, PRL, 89, 134102, (2002)] in case 
of discrete sine-Gordon lattice. 
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Nonlinear phenomena in large diversity of physical sys- 
tems have a close relation with each other because the 
nonlinear dynamics can be described only within few the- 
oretical models Q . Thus there exists a possibility to pre- 
dict novel effects in realistic physical systems via mod- 
eling similar processes in simple hypothetical systems, 
namely chains of coupled nonlinear oscillators could be 
used for this purpose. For instance, as recently has been 
discovered by Geniet and Leon nonlinear supratrans- 
mission phenomenon takes place in discrete sine-Gordon 
lattice, this means that driving harmonically and contin- 
uously one end of the lattice with frequencies within a 
band gap, there is no energy flow through the lattice for 
low amplitude driving, while above definite driving am- 
plitude threshold a sudden energy flow takes place. This 
nontrivial effect has been explained by means of the di- 
rect soliton creation at the end of the lattice, in other 
words the sudden energy flow occurs when the driving 
adjusts the internal oscillations of the localized object. 
It was also noted there the possibility of the existence of 
similar mechanism of gap soliton generation in photonic 
band gap materials. It should be especially mentioned 
that nonlinear supratransmission has been detected not 
only making numerical simulations for model system of 
discrete sine-Gordon lattice, but it is also experimentally 
realized on a mechanical pendulums chain driven at one 
end at band gap frequencies 0- 

The present letter aims to analyze whether the sim- 
ilar scenario takes place in case of discrete nonlinear 
Schrodinger (DNLS) equation and then make the predic- 
tions concerning the corresponding nonlinear processes 
in coupled optical waveguide arrays [j| . The experimen- 
tal conditions are suggested for which optical waveguide 
array becomes transparent with respect to the beam in- 
jected into the single boundary waveguide if beam's in- 
tensity exceeds certain threshold (see FigQJ. This effect 
is opposite to the ordinary case when for low intensities 
(linear regime) the light injected into the single waveg- 



FIG. 1: Suggested experimental setup: The beam is injected 
into the boundary waveguide numbered as j = 0, Qo is a 
linear coupling between boundary and first waveguides and 
Q is a coupling constant between the waveguides in the array, 
no and n are linear refractive indexes of boundary and array 
waveguides, respectively; z is a longitudinal space dimension 
playing a role of time in boundary driven DNLS equation Q . 



uide spreads to other waveguides, while in nonlinear case 
the light is trapped into several neighboring waveguides 
leading thus to the spatial discrete optical breather cre- 
ation (see e.g. most recent experimental papers on the 
subject [4[). It should be just mentioned that longitu- 
dinal space dimension in optical waveguide array plays 
a role of a time variable and one should especially care 
about this suggesting real experiments on waveguide ar- 
rays. 

Let us start form the consideration of boundary driven 
DNLS equation which could be written in the following 
form (j=l, ... , N): 



^+^i+^-i+2|^| 2 ^ = 0; 



■00 = Ae 



■Az 



(1) 



Here z variable stands for the time, A and A are driving 
frequency and amplitude of the boundary. Initial condi- 




FIG. 2: (Color online) Three dimensional plot of time evo- 
lution of boundary driven DNLS equation for inband driv- 
ing frequency A = 1.94 and very small driving amplitude 
Ao = 0.01. As seen the intensity transmits to remote sites. 

tion reads as tpj (0) = and it is supposed that the driving 
is turned on adiabatically, e.g. A = A [l — exp(— z/r)l in 
order to avoid the appearance of the perturbations from 
the initial shock. In simulations it is taken t — 10 and 
nonlinear dynamics is monitored up to the time scales 
10 4 , thus in stationary regime (z 3> r) one has A = Aq. 
The damping also has been applied at the right end of 
the waveguide array in order to suppress edge reflection. 
Note that in the absence of driving the sum of intensities 
Yl j iV'il 2 is a conserved quantity, thus the nonlinear dy- 
namics across the array could be described via intensity 
flux Jj through the site j: 

J,-=z(^* +1 -V*^+i) (2) 

where the intensity and intensity flux at site j are con- 
nected with each other via discrete continuity condition 
d\ipj\ 2 /dz+ (Jj -Jj-i) =0. 

The numerical simulations have been performed choos- 
ing different values of boundary driving parameters A 
and A. From the numerical experiments it follows that 
boundary driving leads to the perturbation of all sites if 
driving frequency A is located within the linear phonon 
band —2 < A < 2, i.e. there is a nonzero intensity flux 
for any driving amplitudes (see Fig. [2J. On the other 
hand if the driving frequency is in upper band gap, par- 
ticularly A > 2, for low driving amplitudes only several 
neighboring sites are excited and intensity flux to remote 
sites is zero. The energy starts to flow only if the driving 
amplitude exceeds certain threshold (see Figs. 0and0}. 
Note that if the driving frequencies are within a lower 
ban gap there is no intensity flow for any driving ampli- 
tudes. Now it is time to analyze the mechanisms for that 
effect. 

As far as a boundary driving is applied it is natu- 
ral to expect that localized solutions will excite. This 
statement is in full accordance with the consideration 
of similar process in discrete sine-Gordon type models 
where the nonlinear supratransmission has been discov- 




FIG. 3: (Color online) The same as in Fig. [5] for band gap 
driving A = 2.04. For driving amplitudes below the threshold 
(A t h — 0.202) the pattern in this graph could be described 
by standing breather solution Q. Inset shows dependence 
of driving threshold upon the driving frequency, asterisks are 
results of numerical simulations on boundary driven DNLS 
equation Q and solid line represents analytical curve ©. 
Dashed line divides the range of Ao (left hand side) for which 
analytical approximate semi-discrete approach given by for- 
mulas and |0. Moreover within the same range above 
threshold low amplitude semi-discrete envelope solitons par- 
ticipate in band gap transmission (see Fig. 0J. While at the 
right side of the dashed line in the inset high amplitude mov- 
ing breathers are excited which are further trapped by the 
lattice (see Fig. |5J, which leads to the suppression of band 
gap transmission. 

ered Q • Thus one can assume that nonzero intensity flux 
will appear when boundary driving excites moving local- 
ized solutions. It is easy to derive a semi-discrete ap- 
proximate envelope soliton solution substituting ansatz 
tpj = ex^>{i(f3z — xj)} into the DNLS equation QJ. 
Then assuming that the envelope ^(j) varies smoothly 
along the lattice and expanding 9(J ± 1) about the site 
j we get the following approximate one soliton solution 
(see for the details of similar derivation e.g. in Ref. |5(): 



, = \p\max i(0z-xj) (3) 

cosh|_|^j| maa ;(j - Vz)\ 

with a following nonlinear dispersion relation for the car- 
rier wave of the envelope soliton (\ varies from to 27r) 

/3 = 2cosx+|V'ilmo*. ( 4 ) 

and V — d/3/dx = — 2sin% is a soliton's group velocity. 
Note that the assumption that soliton envelope varies 
smoothly along the lattice puts the following restriction 
on the soliton amplitude IV'jlmaa; ^ !• 

It is expected that intensity flux appears in the system 
if driving adjusts the nonlinear dispersion relation Q, i.e. 
(3 = A and IV^Imaa; = Ao. In other words, flux is nonzero 
only if one can find such \ that the following condition is 
fulfilled A = 2 cos \ + • Therefore for inband driving 
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FIG. 4: (Color online) The transmission process for the driv- 
ing amplitudes above the threshold, which could be described 
by a train of gap solitons 0. 



FIG. 5: (Color online) Trapping of large amplitude moving 
gap soliton. For relatively large driving amplitudes the soliton 
starts to move but further the lattice traps it. This itself stops 
the transmission process. 



— 2 < A < 2 the nonzero flux appears even for very low 
driving amplitudes Aq, while in upper band gap 



A > 2 



there exists certain amplitude threshold 



A th = VA-2 



(5) 



(6) 



below which (Aq < A t h) there is no intensity flux into 
the system. Instead, only several sites are excited and 
that pattern could be described by a static breather 
solution of Eq. This solution could be derived from 
the general one © for zero velocity V — — 2 sin \ = 0, i.e. 
X = and requiring the fulfillment of boundary condition 
ipa = Aq exp(iAz) in stationary state (z 3> r): 



VA - 2 



cosh[(j + xo)\/A-2j 



(7) 



where xo — acosh[VA — 2/Ao]/\/A — 2. 

The results of numerical simulations are fully explained 
by the above consideration. In inset of Fig. [31 it is pre- 
sented the comparison between numerical experiments on 
Eq. and analytical formula © for the driving thresh- 
old above which (Ao > Ath) a nonzero intensity flux ap- 
pears in the system. For driving frequencies A close to 2 
there is a perfect agreement, but this agreement becomes 
worse for larger driving frequencies. The point is that 
for driving frequencies sufficiently larger than 2 thresh- 
old amplitudes become comparable with unity according 
to the relation JJJJ. But for such amplitudes the contin- 
uum envelope approximation J3J is invalid, moreover, as 
numerical simulations show large amplitude excitations 
are trapped by the lattice, (see Fig. 0), as a result lo- 
calizations do not move and intensity flux becomes zero. 
Thus the band gap transmission effect exists if there are 
moving solutions in the system. As a result in case of 
DNLS equation the discovered phenomenon is observable 
for driving frequencies 2 < A < 2.09. 



Now let us discuss how the obtained results could be 
applied to describe nonlinear transmission processes in 
the system of coupled optical waveguides. For the pur- 
pose to realize a band gap driving it is suggested (see Fig. 
^) to inject a beam into boundary waveguide with linear 
refractive index no larger than refractive index n of other 
waveguides forming the array. Let us introduce a linear 
coupling constant between array waveguides as Q, while 
the coupling between boundary (j = 0) and first waveg- 
uides is defined as Qo- Besides that, let us suppose for 
simplicity that the nonlinear refractive index (kerr non- 
linearity) in j = waveguide is equal to zero and onsite 
nonlinear refractive index in array waveguides is D. Thus 
the wave envelopes in waveguides could be described by 
a set of the following equations (j = 2, . . . , N): 



«^ 1- — n ipo + QoWi = 

Oz c 



«? 1 + -^i+QoV'o + W2 + ^|V'i| 2 V'i =0 (8) 
oz c 

.dp_ + u + + ^ j + £>|-0j| 2 Vj = 0. 

oz c 

where w is a carrier wave frequency and c is a light 
velocity. Last equation from the set JHJ is a well known 
representation of infinite waveguide array by the DNLS 
equation 0] while first two equations describe influence 
of boundary waveguide (with different linear refraction 
index) on the semi infinite array (see also Ref. 0). After 
the appropriate rescaling 



^=^e iz " n l c ^2QlD for j = l, 



(9) 



tAo = V e* Wc (Q/Qo) ViQ/D; z = z'/Q 

Eqs. JSJl obtain simpler form (A = — nj/Qc): 



dib' 

^ + (^+i + ^-i) + K-l 2 ^ = o 
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FIG. 6: (Color online) Rescaled waveguide intensity \ip'j\ 2 in 
the band gap transmission regime. Inset shows longitudinal 
dimension z' dependence of rescaled beam intensity in bound- 
ary waveguide. Creating the solitons the intensity from the 
boundary waveguide is transferred to other waveguides and as 
a result the transmission stops when beam intensity in bound- 
ary waveguide goes below the threshold. The simulations are 
done for Qo/Q = 0.1. 

(j = 1,...,N) which reduces to the boundary driven 
DNLS ||TJ with boundary condition tp' Q = i/j' (0) exp(iAz') 
in the limit (Qo/Q) — > and therefore for Qo/Q <C 1 
one can use the results derived for the case of boundary 
driven DNLS JSJ and ©. Particularly, for 

A = uj(n - n)/Qc > 2 (11) 

the localized excitations © or 0) form with propagation 
constant located in the upper band gap (3 — A > 2. 
Thus if the injected intensity in the boundary is below 
the threshold 

lM ° )lth -^J^ mth ^DQl[ c 2Q 

(12) 

one has a static breather solution (JJJ) and intensity flux 
to the array waveguides is equal to zero, while above the 
threshold energy transmission begins via gap solitons J3Jl. 

The above expression for the threshold 1|12|) becomes 
exact in the limit (Qo/Q) —> 0. For small but nonzero 
(Qo/Q) the beam intensity in the boundary waveguide 
could be considered as almost constant quantity irrespec- 
tive to the spread of energy in nonlinear band gap trans- 
mission regime, because according to the rescaling beam 
intensity in the boundary waveguide is Q/Qo times larger 
than the amplitude of propagating soliton through the ar- 
ray. In case when Qo becomes comparable with Q almost 
all intensity in boundary waveguide is needed to form a 
soliton. As a result the intensity in boundary waveguide 
sharply decreases and therefore much larger threshold in- 
tensity is needed [than given by relation l)12[l] to develop 
the nonlinear transmission in band gap regime. In Fig. 



|B] it is presented the picture for the case Qo/Q = 0.1. 
As seen the beam intensity in boundary waveguide is 
above the threshold at the origin z' = and it produces 
gap solitons causing the band gap transmission. But this 
process itself causes the decrease of beam intensity in the 
boundary waveguide, the intensity after the creation of 
several solitons goes below the threshold and transmis- 
sion process is not observable for large z' . 

As it has been mentioned above nonlinear band gap 
transmission regime holds if one has low amplitude soli- 
tons. Large amplitude solitons tend to pin and energy 
transfer becomes much less effective. However, this is 
true only in case of DNLS equation where only onsite 
nonlinearities are taken into account. As shown recently 
[8j considering also the terms describing also intersite 
nonlinearities one has moving breather solutions even 
at large excitations amplitudes. High intensity moving 
breathers have been also detected on the recent experi- 
ments Q . Therefore the numerical simulations have been 
undertaken adding to DNLS equation also terms with in- 
tersite nonlinearities. In this case the energy transfer via 
the moving breathers take place even for large excitation 
amplitudes and as a result optical transparency regime 
is observable in whole range of A > 2. 

Summarizing it should be noted again that the novel 
scenario of nonlinear band gap transmission in optical 
waveguide arrays is predicted and simple experimental 
setup is suggested for its realization. Suggested experi- 
mental setup would serve also for generation of optical 
gap solitons propagating across the waveguide array. 
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